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1. Introduction 



The confinement puzzle has been with us ever since the birth of quantum chromodynamics 
(QCD). By means of lattice calculations, it has been possible to penetrate the infrared nonpertur- 
bative sector of QCD and recover a confining potential between (static) quarks At present, 
however, available lattice sizes do not suffice to describe the Green functions in the deep infrared 

The continuum approach, on the other hand, has the intriguing feature that the infrared limit 
can be studied asymptotically. In the last decade, a new understanding of infrared QCD has arisen 
from studying continuum Yang-Mills (YM) theory via Dyson-Schwinger equations. The Landau 
gauge has the advantage of being covariant and therefore encouraged many to intensive investiga- 
tion of the infrared properties of YM theory [Q,^. In Coulomb gauge, non-covariance brings about 
severe technical difficulties which are only recently on the verge of being overcome [^. Neverthe- 
less, the Coulomb gauge might be the more efficient choice to identify the nonabelian degrees of 
freedom. It is well-known that screening and anti-screening contributions to the interquaik poten- 
tial are neatly separated in Coulomb gauge perturbation theory As for the infrared domain, the 
Gribov-Zwanziger scenario serves as a transparent confinement mechanism ||]. 

A further advantage of working in the physical Coulomb gauge is that one may pass over 
to a Hamiltonian description. This opens up direct access to the heavy quark potential via the 
expectation value of the Hamiltonian. In recent years, variational methods have been pursued to 



solve the Yang-Mills Schrodinger equation with a Gaussian type of wave functional M, 10, 11 



12]. Despite Feynman's critique [13], it turns out that the wave functional is sensitive to infrared 



modes and the variational method a powerful tool, at least for the qualitative description of YM 



theory. With careful treatment of the operator ordering in Coulomb gauge [|14|], it is possible to 
find a strictly linearly rising heavy quark potential. We report on the latest results found in the 
Hamiltonian approach to YM theory in Coulomb gauge. This includes the full calculation of gluon 
and ghost Green functions and a running coupling. Furthermore, the 't Hooft loop, an (dis-)order 
parameter for confinement, will be calculated using the results of the Green functions. 

This paper is organized as follows. In section ^, the Yang-Mills Hamiltonian in Coulomb 
gauge and the equations of motion are introduced. The latter will be solved variationally and the 
heavy quark potential and the running coupling are presented in section |[ The 't Hooft loop is 
discussed in section ^ and conclusions are given in section ^. 



2. Yang-Mills Schrodinger equation and Dyson-Schwinger equations 

In the canonical quantization approach, we choose Aq{x) = and impose the usual equal-time 
commutation relations among the gauge field A"{x) and the conjugate momentum Tl"{x) to arrive 
at the Weyl gauge Hamiltonian. Since Aq originally serves as the Lagrange parameter of the Gauss 
law, the choice of Weyl gauge requires a restriction on the Hilbert space, 

Dn\^')=gp,„\^') (2.1) 

where g is the gauge coupling, Pm(x) the density of external color charges, and D"^ = diS"^ + gA"}^ 
with A?^ = A'jf'^''. Fixing the residual time-independent local gauge invariance by the Coulomb 
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gauge, dfAi = 0, and eliminating the longitudinal part of the momentum operator n by means of 



the Gauss law ([2.1|), one arrives at the Hamiltonian that depends only on transversal fields, 

H = ^J{^-'n^n + B^ + g^^-'pFJ^p). (2.2) 

The appearance of the Faddeev-Popov determinant ^[A] = det{—D[A]d) is due to a non-trivial 
change of coordinates to the transverse fields and turns out to be crucial to the infrared properties 



of the theory. The latter term in Eq. ( [2.2[ ) describes the Coulomb interaction of dynamical and 
external charges p = —An + p,„ via 



F"'{x,y) = {x,a\ {-Ddy' {-d') {-Ddy' \y,b) (2.3) 

and reduces to the familiar Coulomb law in the abelian theory. 

With the Hamiltonian (p!2|) at hand, we may apply the variational principle to find the wave 



functional *P[A] = (A|*F). Inspired by QED, we choose [10] 



^[A] = / ^Aexp ( - / A(oA ] (2.4) 



with a normalization constant The factor of ^ is chosen to alleviate the computation of 
expectation values, similar to defining radial states in quantum mechanics. A different power of 



^ in the wave functional does not change the properties of the solution [15]. One may think of 
the variational parameter CO as in the inverse of the gluon propagator, 

Dlfix,y) = Af(^)A^(j) \^) = ^-5"%ix)co-\x,y) (2.5) 

with tij being the transverse projector. It is determined by solving the functional Schrodinger 
equation, i.e. minimizing the energy H This gives rise to a non-linear integral equation in 



CO which we refer to as the gap equation. It was derived to two-loop order in the energy in Ref. [jlOp 
and reads in momentum space {k= [k[) 

0)^{k)=k^ + xHk)+Uk)+I^a,- (2-6) 
Here, x (k) abbreviates the so-called curvature and it is related by 



to the ghost propagator 



{^\ i-Dd)-' = . (2.8) 



A Dyson-Schwinger equation for the ghost form factor d may be derived from the path integral, or 
alternatively from the following operator identity for G[A] = (— D5)^\ 

G[A] = {-BY' + {-d'^y'gAd G[A] (2.9) 
d-Hk)-,-'-^ f (I (k g)^) ^(l'^-ql) (210) 



which yields 
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Both in the ghost Dyson-Schwinger equation ( 2.10| ) and in the equation for the curvature (2.7), we 



have approximated the proper ghost-gluon vertex by its tree-level counterpart Y^. This amounts to 
the factorization 

{AiG[A]) = Dij {G[A] r']G[A]) « A7 {G[A])r^j {G[A]) . (2.11) 
The non-renormalization of the ghost-gluon vertex in gauges where the gluon propagator is trans- 



verse, such as the Coulomb and the Landau gauge [[16|, |17|], suggests that the above approximation 
is a good one. Dyson-Schwinger studies in both four and three-dimensional Landau gauge as well 
as lattice calculations in four-dimensional Landau gauge confirmed that the dressed vertex is close 
to tree-level [Il8]|. The case of three-dimensional Landau gauge resembles the Coulomb gauge and 



therefore we adopt the approximation (2.11). This vertex' non-renormalization will have crucial 
impact on the IR sector of the solutions. 



The other momentum dependent term 1 01 (k) in the gap equation ( |2.6| ) reads 
f (,._^c f d^q 2^ d{k - q) V(k - q) [(o{q) -z{q) + Z(k)]' - (o{kf 

and is due to the Coulomb interaction part of the Hamiltonian. Here, the form factor / measures 
the deviation from the factorization of the Coulomb potential, 

{G[A]{-d')G[A]) = {G[A]) {-d')f {G[A]) . (2.13) 

In the infrared, we set f{k) = 1, factorizing the expectation value for the Coulomb propagator 



(2.13) equivalently to the one for the ghost-gluon vertex in Eq. (2.11). In the ulti"aviolet, f{k) is 



treated perturbatively, see ||10|]. 

In order to fix the Coulomb gauge uniquely, configuration space must be restricted to the 
compact fundamental modular region. As suggested in [^], this entails the "horizon condition" for 
the ghost form factor, 

d-\0)=0. (2.14) 

As we shall see, the horizon condition ( 2.14[ ) has the consequence that all form factors d, x and (O 
diverge in the infrared. 

3. Green functions, heavy quark potential and running coupling 

The ultraviolet divergences encountered in the gap equation ( |2.6| ) are removed by subtracting 
the equations at an arbitrary renormalization scale jj.. Alternatively, one can eliminate the diver- 



gences by adding appropriate counter terms to the YM Hamiltonian and to In ^ [|19|]. This elim- 
inates the UV-divergent constant /J^ from Eq. ( ^^ and involves some renormalization constants, 
one of them can be chosen as c = lim<:^o(w(^) — and fixed by the requirement of minimal 
energy to be c = 0. For details, see Ref. [p4|]. 

The solutions for the form factors 0){k), d(k) and xiki) can be seen in Fig. [!} In the asymptotic 
infrared, the gluon form factor ft)(^) approaches the curvature xik), reflecting the dominance of 
ghost degrees of freedom, cf. Landau gauge [Q]. With x{)i) being infrared enhanced, the ghost 
content of the solution makes propagation of gluons over asymptotically large distances impossible, 
hence gluons are confined. 
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Figure 1: Left: the gluon form factor co{k) and the curvature Right: the ghost form factor d{k). 



The correlation of the asymptotic infrared power laws is due to the non-renormalization of the 
ghost-gluon vertex, 



(3.1) 



Without imposing the horizon condition ( |2.14| ), solutions to d, X and CO can be found that approach 
finite values in the infrared [ |l9| ] so that the energy functional is dominated by the ultraviolet modes, 
as speculated by Feynman [pj|]. Conversely, the infrared power law solution ( |3.1| ) where the horizon 
condition is satisfied are not subdominant to ultraviolet modes and turn out independently of the 
details of the wave functional. Even a stochastic vacuum, *P[A] = 1, would produce the same results 
for the infrared [|l^]. One may therefore be confident using the variational principle. 

The infrared enhancement of the form factors co{k) and d{k) is qualitatively reproduced by 
recent lattice calculations [^]. 

Equipped with the ghost form factor d{k), the heavy quark potential can be found by choosing 

p^(x) = 5"' (5^'\x-r/2) - 5'^'Hx + r/2). 



(3.2) 

and recalculating the energy {H) with fixed CO. There is only one contribution to the energy that 
depends on the distance r between the quarks. Using Eq. ( 2.13| ), it reads 



d'q d\q)f{q) 
(271)3 ^2 



;i -.'■"■■•). 



(3.3) 



With the infrared behavior of the form factors ( |3.l| ), we find that Vc{r) rises linearly in the infrared 
and thus confines heavy quarks as shown in Fig. ^. By matching the slope of the hnear potential to 
the lattice string tension a, one may set the scale. 

Apart from the solution in Fig. |l| with the asymptotic infrared behavior ( \A ) there is one further 
solution with slightly different infrared exponents for the power laws. The latter was discovered 
first [10], however, it does not have the same attractive features as the one in Fig. |TJ In particular, 
the heavy quark potential is strictly linearly rising only for the solution presented here. 

A nonperturbative running coupling may be extracted from the ghost-gluon vertex |[T^, In]], 



aik) = —kd^(k)(o-\k), 
3n 



(3.4) 
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Figure 2: Right: the Coulomb potential V{r). Right: the running coupling a{k). 
With a tree-level vertex, it can be shown that one finds a finite value in the infrared, a(0) 



16;r/(3A'c) [11 1. In the ultraviolet, we find the correct \/ln{k/[j.) scaling from one-loop pertur- 
bation theory. However, the first coefficient of the beta function, j3o, is off by a factor of 8/11. 
This is due to the approximations and requires further investigation. Recent lattice calculations 



in Coulomb gauge [21] show qualitative agreement with the running coupling presented here. In 
comparison to the analogous running coupling in the Landau gauge [Q], we do not find a bump for 
intermediate momenta yielding a spurious zero in the beta function. Note that a{k) in Fig. § is a 
monotonic function. 

4. The 't Hooft loop 

A (dis-)order parameter of confinement is the 't Hooft loop (V(C)) [ p2| ] whose operator V(C) 
is defined by the relation V{Ci)W{C2) = Z^'^'^'-'^^'^W{C2)V{Ci), where W{C) is the operator of the 
spatial Wilson loop, Z is a (non-trivial) center element of the gauge group and L(Ci,C2) denotes 
the Gaussian linking number. An explicit realization of V{C) in continuum Yang -Mills theory was 



derived in ref. [23] and is given by 

V(C)=exp i f d^x£/i"[C]ix)n'^{x) 



(4.1) 

Here £/[C] denotes the gauge potential of a (spatial) center vortex whose magnetic flux is localized 
at the loop C. Since V{C)^{A) = *F(A + i/[C]) the 't Hooft loop is a center vortex generator. 
Using the wave functional found in the variational solution of the Yang-Mills Schrodinger equation 
in Coulomb gauge, as described above, the expectation value {V{C)) = exp(— ^(C)) was evaluated 
for a planar circular loop C and is was found that the exponential S{C) obeys a perimeter law 



signaling confinement []24[]. This result is in accord with the linear behavior found for the static 
color potential. 

5. Conclusions 

We have solved the Yang-Mills Schrodinger equation approximately and thus determined the 
vacuum wave functional. Our solutions exhibit the phenomena of confinement of gluons as well 
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as heavy quarks. As an improvement on previous results, the heavy quark potential rises strictly 
Unearly. The nonperturbative running coupling derived from the ghost-gluon vertex was presented 
and the 't Hooft loop was calculated. It is promising that the results have the crucial features of 
nonperturbative physics, and that calls for further investigations in the Hamiltonian approach. 
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